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We present an ah initio study of the effect of macroscopic deformation on energetics of twelve 
alloying elements in bcc Fe under three specially designed strain modes. We find that there exists 
a universal linear relation of describing the volume dependence of substitutional energy of alloying 
elements via introducing two factors — the system shape factor (/ss) and the size factor of alloying 
element M (fig/): E S ub ~ /ss^g/ V. fi^ well describes the effect of intrinsic alloying-atom size and 
the influence of chemical interaction with matrix atom, and / ss characterizes the degree of system 
lattice distortion under deformation. This relation is further validated using the published data of 
stained-modulated doping in GaP [Phys. Rev. Lett. 105, 195503(2010)]. 
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Fe-bascd alloys are among the most widely used ma- 
terials, of which the ferritic martensitic steels represent 
a technologically important class of materials with many 
applications in fission and fusion energy systems. Alloy- 
ing solutes in Fe-based alloys are of great importance in 
modifying physical properties, and particularly their con- 
tents strongly govern the mechanical performances such 
as the resistance to hardening and embrittlement under 
neutron and proton irradiation [ll-Q- It is found that 
the optimal Cr content to improve one special property 
may be different with that to improve another for bcc Fe 
based alloys, e.g., the content of 2 to 6% Cr to reduce 
the irradiation swelling whereas 9% Cr to reduce ductile 
brittle transition temperature (sec Ref. 3 and references 
therein). Most previous investigations of alloying ele- 
ments in ferritic martensitic steels have been performed 
in macroscopically undeformed crystals — free of any 
macroscopic stain, except those induced by the solutes 
or defects themselves. Nonetheless, as the structural ma- 
terials in fission and fusion reactors, the ferritic marten- 
sitic steels experience both multiaxial loading and fairly 
high temperatures, and the internal strain field becomes 
more complex due to the presence of radiation damages. 
In fact, all materials usually undergo macroscopic de- 
formations due to externally applied loads or in special 
in-service conditions. So in recent years many researches 
have been made on the properties of solutes and point de- 
fects under macroscopic deformations based on electronic 
structure calculations. Gavini found that the volumet- 
ric strain associated with a deformation largely governs 
the formation energies of monovacancies and divacan- 
cies in Al and concluded the nucleation of these defects 
is increasingly favorable under volumetric expansion 
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Chen et al. reported that the strain has remarkable influ- 
ence on the stability, reorientation and migration of self 
interstitial atoms (SIA) in bcc Fe [f|: for instance, uni- 
axial expansion induces a SIA spontaneous reorientation 
from < 111 > to < 100 >. Zhu et al. found in zincblende 
GaP the impurity formation energy changes monotoni- 
cally in a linear fashion with the applied external strain 
and thus proposed that the intriguing unbounded strain- 
induced change in impurity formation energy can be used 
effectively to enhance dopant solubility in a wide range 
of semiconductors Q. These results raise an interest- 
ing question on how macroscopic deformation affects the 
stability of alloying elements in Fe-based alloys, which to 
our knowledge remains open. If the macroscopic defor- 
mation does influence on the energetics of alloying ele- 
ments, the aggregation and precipitation may take place 
consequently, and thus the mechanical performance may 
be degraded. 

In this Letter, we present an ah initio study of the ef- 
fect of macroscopic deformation on substitutional energy 
of twelve alloying elements (Al, Co, Cr, Cu, Mo, Nb, Ni, 
Si, Ta, Ti, V, and W Q) in bcc Fe under three different 
strain modes. Hydrostatic strain mode is to explore the 
pure volume change effect, and a set of submodes in both 
normal strain and shear stain modes are to explore the 
lattice distortion effect. We find that the macroscopic de- 
formation effect on substitutional energy of alloying ele- 
ments is governed linearly and universally by two factors, 
the system shape (lattice distortion) factor and the size 
factor of alloying element. The size factor of alloying el- 
ement well describes the effect of intrinsic alloying-atom 
size and the influence of chemical interaction with ma- 
trix atom, and the system shape factor characterizes the 
degree of lattice distortion under deformation. 

The present calculations are performed within spin- 
polarized density functional theory as implemented in the 
Vienna ah initio simulation package (VASP) Q. The in- 
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TABLE I: Strain components (%) of each submode in NS and 
SS. 

en, e 22 , e 33 (NSl-4)/eu, e 33 , e 12 (SS89-85) 



Vol. contracted 
by 3.0% 



Vol. unchanged 



Vol. expanded 
by 3.0% 



NS1 

NS2 

NS3 

NS4 

SS89 

SS88 

SS87 

SS85 



-1.98, - 
-2.94, 
-3.88, 
-4.81, 
-0.996, 
-0.985, 
-0.966, 
-0.906, 



-1.0, - 
-1.0, 
-1.0, 
-1.0, 
-1.0, 
-1.0, 
-1.0, 
-1.0. 



-0.01 
0.98 
1.97 
2.96 
0.864 
1.73 
2.59 
4.33 



-0.99, 
-1.96, 
-2.91, 
-3.85, 



1.0 
2.0 
3.0 
4.0 



0.004,0, 0.873 
0.015, 0, 1.75 
0.034, 0, 2.62 
0.096, 0, 4.37 



0, 1.0, 2.01 
-0.98, 1.0, 3.02 
-1.94, 1.0, 4.03 
-2.88, 1.0, 5.04 
1.004,1.0, 0.881 

1.02, 1.0, 1.76 

1.03, 1.0, 2.65 
1.10, 1.0, 4.41 



teraction between ions and electrons is described by the 
projector augmented wave (PAW) method @. Exchange 
and correlation functions are taken in a form proposed 
by Perdew and Wang (PW91) [HJ within the generalized 
gradient approximation (GGA). The supercell approach 
with periodic boundary conditions is used. The super- 
cell contains 127 Fe atoms and one alloying atom. We 
used an energy cutoff of 350 eV for the plane-wave ex- 
pansion of the wave functions and a 3 x 3 x 3 fc-point mesh 
for Brillouin zone sampling according to the Monkhorst- 
Pack scheme. All atomic relaxation calculations are per- 
formed at constant volume and shape using the conjugate 
gradient algorithm. For simplicity the side length of the 
forcmentioned bec cubic supercell in equilibrium is set to 
be four times of the optimized lattice constant of bec Fe 
(2.833 A). In the calculation of substitutional energy the 
following formula is used: 



E su b = E(nFe + 1M) 



n + 1 



xE((n+l)Fe)-E(M), (1) 



where E(nFe + 1M) is the total energy of a bec cubic 
supercell containing n Fe atoms and one alloying atom 
M , E((n + l)Fe) is the total energy of the same bec 
cubic supercell filled with only Fe atoms under the same 
deformation, and E(M) is the energy per atom of pure 
crystal M in its stablest phase. 

In this work, the macroscopic deformation is described 
by a macroscopic strain tensor with the coordinate axis 
1, 2 and 3 being chosen along [100], [010] and [001] di- 
rections of the perfect bec crystal, respectively. Because 
of strain tensor belonging to a six-dimensional space, a 
complete characterization of deformation effect on the 
energetics of alloying atoms is beyond reach. Here we fo- 
cus on three different modes of strain: hydrostatic strain 
(HS), normal strain (NS) and shear strain (SS) within 
relative volume variation from —3.0% to 3.0%. Both 
NS and SS modes include four submodes: NS1-4 and 
SS89-85, as listed in Table [J Here, for example, SS89 
denotes that the angle between two adjoining surfaces is 
changed to 89° from initial value 90°. Note that, 1) in 
NS/SS mode the cube symmetry of the system changes 
into orthogonality/monoclinism; 2) in each submode of 



NS and SS a set of strains arc applied in such a way 
that leads to different volumetric change but the sys- 
tem being kept similar shape, in other words, the ra- 
tio of three edges and the angle between any two inter- 
sected edges of the supercell are unchanged. For exam- 
ple, in submode NS1 the system is in rectangular paral- 
lelepiped shape and the ratio of three different edges is 
fixed at 0.99 : 1.0 : 1.01. From NS4 to NS1 (from SS85 
to SS89) the system shape approaches the perfect cube 
step by step. For all submodes of SS mode the follow- 
ing equations exist: en = e22, e23 = ^31 = 0. Of these 
submodes the strain components with relative volume 
changes of — 3.0%(contracted), 0%(kept at equilibrium) 
and 3.0% (expanded) are summarized in Table HI 
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FIG. 1: (Color online) Substitutional energy versus volume 
for twelve different alloying elements in HS mode. Symbols 
represent the calculated results and solid lines are the corre- 
sponding linear fits. 

We first investigate the isotropic volumetric deforma- 
tion influence on substitutional energy of alloying ele- 
ments in HS mode. Figure[T] displays the substitutional 
energy as a function of volume ranging from 0.97Vo to 
1.03Vo in HS mode. In our present work Vq is the equi- 
librium volume of the supercell with 128 Fe atom. As 
can be seen, the substitutional energy of Si is insensitive 
to isotropic deformation, while with the increase of vol- 
ume those of Al, Co, Cr, Cu, Ni, and V decrease slowly 
and those of Mo, Nb, Ta, Ti and W decrease strongly. 
The results suggest that in HS mode the stability of the 
alloying elements studied except Si are considerably influ- 
enced by hydrostatic volumetric deformation, and these 
alloying atoms except Si prefer strongly to stay in the 
expansion zone, which in turn changes the distributions 
of alloying atoms in bec Fe, considering the direct corre- 
lation between substitutional energy and solubility. For 
every studied alloying element the substitutional energy 
depends linearly on the volume with a negative slope 
(except Si with a much small positive slope), which is 
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dependent on alloying element. Interestingly, when ex- 
amining the calculated substitutional energies of each al- 
loying element in each submode of NS/SS with volumes 
of 0.97Vb, Vq, and 1.03Vo, we find the similar linear vol- 
ume dependence as in HS mode but with an increased 
slope. The obtained linearity has been further checked 
via the added calculations for Al, Cr, and W in NS1, 
NS4, SS89, and SS85 submodes with volume of 0.982Vo, 
0.994Vb and 1.018F , and here only the results of W in 
NS4 and SS85 submodes are presented in FIG.[2] Fig- 
ure^] clearly indicates that, there does exist a similar lin- 
ear volume dependence of substitutional energy in sub- 
modes in NS and SS modes as in HS mode but with a 
larger slope. For example, in HS the isotropic volumetric 
increase from 0.97Vo to 1.03Vo leads to a substitutional 
energy decrease of 0.234 eV for Cr and 0.575 eV for W, 
while in NS4/SS85 the substitutional energy decrease are 
0.162/0. 116 eV for Cr and 0.486/0.442eV for W. Note 
that AE su b (= S su b|y=i.o3V — £sub|v=o.97Vo) divided 
by 0.06Vb is approximately identical with the element- 
dependent slope of the forementioned linearity between 
the substitution energy and volume. These results tell us 
that in any submode of NS/SS the volumetric deforma- 
tion effect on the stability of the studied alloying elements 
except Si is noticeably reduced compared with that in HS 
mode. Our observed linear volume dependence of substi- 
tutional energy are in good agreement with the reported 
linear behavior of impurity formation energy in semicon- 
ductors under hydrostatic strain by Zhu et al Q. They 
found that the change of impurity formation energy is a 
linear function of volume and the slope is proportional 
to the effective size difference between the host and the 
dopant atom; and trends of the strain-induced change of 
impurity formation energy for the biaxial and the hydro- 
static strain are about the same, except that the effect 
of the biaxial strain is smaller. Thus, the deformation 
in NS and SS modes can modulate to some extent the 
volume dependence of substitution energy compared to 
pure isotropic volumetric change in HS mode. 
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FIG. 2: (Color online) Substitutional energy -E su b of W versus 
volume in HS, NS4 and SS85. Lines are linear fits to the data. 

Above results reveal for every studied alloying element 
the existence of linear volume dependence of the substi- 
tutional energy in the same deformation mode (i.e., the 
system is kept a similar shape such as in the HS mode 
and in the submode of NS/SS mode), but the slope is 



dependent on the alloying element. As can be seen from 
FIG.[T] for the alloying atom with much larger radius 
the slope is smaller than that with a smaller atomic ra- 
dius (e.g., Ta and Si). This implies that the difference in 
atomic size between the alloying atom and the host is a 
key factor for further exploring the linear volume depen- 
dence of the substitutional energy. Recently, Olsson et 
al. proposed a special method to calculate the size fac- 
tor (fi^r ) of single solute atom M in bec Fe matrix, and 
found that the calculated fi^f can reflect not only the 
effect of the solute atom size but also the influence of the 



chemical interaction with the matrix [ill . Following this 
method, we first obtain the bulk modulus B for every bi- 
nary alloy by fitting the linear relation between pressure 
and volume (A_P = — BAV/Vq) in HS mode with volume 
varying from 0.97Vb to 1.03Vb- Then the size factor of 
the alloying atom is calculated by 



M 



-VoAP ex /(Q FeJ B), 



(2) 



where AP ox is the excess pressure of the binary alloy and 
Q,F e is the atomic volume of bec Fe. AP ex is calculated by 
subtracting the external pressure of pure bec Fe matrix 
in HS from that of the binary alloy in HS with the same 
equilibrium volume Vq. If the 17^/ of alloying atom cor- 
rectly describe the intrinsic size difference and the change 
in atomic interaction, then there may exists a general (or 
universal) law of the macroscopic deformation effect on 
substitution energy. So we replot the substitution energy 
changes AP su b of twelve alloying elements as a function 
of in FIG. [3] To our surprise, as displayed in FIG. [3] 
it is found that in those modes/submodes whose system 
being kept a similar shape AP su b exhibits a universal 
dependence on the Cl^ : for all studied alloying elements 
AP su b (i.e. the slope of linear volume dependence of 
Psub) is approximately inverse proportional to fi^f. In 
Ref. 6, it is obtained as a first order approximation that 
Pdop^a(Vhost-Vhost+do P ant)^, where a is the elastic con- 
stant, Vhost+dopant is the equilibrium volume of the sys- 
tem with the dopant and Vhost is the equilibrium volume 
of the host lattice. This formula can also well describe the 
linear relation of substitutional energy versus volume in 
FIG.Q] and the slope is determined by a(Vh os t — Miost+M)- 
Based on the above discussion a(V"host — ^4ost+dopant) m 
Ref. 6 corresponds to Q^f in present work. Therefore, 
using Eq.[5] (where AP CX is replaced by — BAV ex /Vo and 
f^Fe is replaced by f^Ga) and the data summarized in 
Table I in Ref. 6 on strain-enhanced doping in semi- 
conductors, we calculate f2^ of doping elements in GaP. 
The doping energy difference APdop resulted from the 
strain change between -2% and 2% in HS mode as re- 
ported in FIG. 1 in Ref. 6 is redisplayed as a function 
of the calculated fifj? in the inset of FIG. 3. This in- 
set figure clearly indicates that for doping of zincblende 
GaP the doping energy changes resulted from hydro- 
static strain show a similar universal dependence on the 
O^, though the doping energies of these elements ex- 
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TABLE II: The system shape factor of each NS/SS submode 
in unit of eVA 3 . The system shape factor in HS mode is 
-2.51 eVA" 3 . 



Co V Ni Ti 



Mo WTa Nb 



Mode 


NS4 


NS3 


NS2 


NS1 


SS85 


SS87 


SS88 


SS89 


/ss 


-2.31 


-2.38 


-2.40 


-2.48 


-2.39 


-2.40 


-2.45 


-2.48 



hibit different, even opposite volume dependence. Thus 



the Q, 



M/D 



is a suitable factor which well describes the 



intrinsic size difference and the change in atomic interac- 
tion between alloying/doping atom with host atom and 
produces a universal effect of pure volumetric deforma- 
tion. In addition, as forementioned and shown in FIG. 2 
the element-dependent slope of the linearity between the 
substitution energy and volume is larger in any submode 
of NS/SS than in HS mode, suggesting that the system 
shape factor (/ ss ) should be introduced to describe the 
lattice distortion of system under strained conditions. So 
the substitutional energy change due to macroscopic de- 
formation can be described universally by the following 
formula: 



Ungv + 1. 



(3) 



Here, the constant /intercept / is dependent on both the 
alloying element and the strain submode. By linear fit- 
ting of the calculated Ai? su b versus the / ss of each 
NS and SS submode is obtained and listed in Table HI1 As 
can be seen, from NS4 to NS1 and from SS85 to SS89 the 
system lattice distortion is weakened, the system shape 
factor values become smaller ant approaches the value in 
HS mode where no lattice distortion. This result suggests 
that the system shape factor can well describe the system 
lattice distortion. Therefore, the macroscopic deforma- 
tion effect on substitutional energy can be divided into 
two parts: the system shape distortion effect described 
by the / ss factor and the pure volumetric change effect 
generalized by the for different alloying elements. 
In shape/lattice distorted systems the pure volumetric 
change effects on substitutional energy is enlarged for 
the alloying element who has negative size factor, and is 
reduced for the one who has positive size factor compared 
with the cases in perfect cubic systems. 

In summary, the ab initio calculations reveal that 
within volume change ranging from 0.97Vo to 1.03Vo 
the macroscopic deformation effect on the substitutional 
energy of alloying elements can be divided into pure 
volume change effect and lattice distortion effect and 
be well described by a universal and simple relation: 
E su b ~ fmQtfV. This relation should be beneficial to 
understanding the mechanical performance degrading of 
Fe-based alloys under deformation. We believe that this 
relation can be generalized to other systems in studying 
substitutional or doping energy change under deformed 
condition. In addition, in the lattice distorted systems 
the pure volume change effects on substitutional energy 
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FIG. 3: (Color online) Substitutional energy change Ai5 su b 
versus size factor in HS, NS4 and SS85 with volume change 
from 0.97Vo to 1.03 Vb for twelve alloying elements. Inset: 
The doping energy difference resulted from the strain changes 
between —2% and 2% in HS mode in GaP versus size factor 
of dopant. Lines are linear fits to the data. 



is enhanced for those alloying elements who have negative 
size factor, and is reduced for those who have positive size 
factor compared to the cases without lattice distortion, 
which may be instructive to design better materials. 
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